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Abstract 


The refilling of propellant tanks while in a low-gravity environaent requires that en- 
trapped vapor bubbles be collapsed by increasing the systea pressure. Tests were per- 
foraed to verify the aechanism of collapse for these large vapor bubbles with the therao- 
dynaaic conditions, geoaetry, and boundary conditions being those applicable to propellant 
storage systems. For these conditions it was found that conduction heat transfer deter- 
mined the co’ apse rate, with the specific bubble geoaetry having a significant influence. 


Introduction 


The capability of refilling the storage tanks for liquid propulsion systems while a 
spacecraft is in earth orbit will yield a number of benefits. With refillino, the life of 
a spacecraft can be extended or a space-based vehicle could be reused for rumerous mis- 
sions. Methods of refilling a pr-ipellant tank while in low-qravity are currently being 
developed. This paper considers one aspect of tank refilling: the collapse of vapor bub- 

bles that may become entrapped within the tank during filling. 


Expulsion of liquid propellants from a tank under low-gravity conditions requires some 
means of ensuring that only liquid will be supplied to the engine. Capillary propellant 
aanageaent devices, using fine-aesh screen to orient liquid and exclude gas are the likely 
choice for the expulsion system, especially for cryogenic propellant applications. These 
devices are now being sed for propellant expulsion on the Space Shuttle^ and various 
communications satellites. In one configuration, the fine-mesh screen is used on channels 
moui ad near the tank wall and encircling the tank. These channels form a flow passage 
from the bulk liquid, regardless of its orientation, to the tank outlet. Liquid flows 
through the screen in preference to gas, due to the capillary pressure differential devel- 
oped at the pores of the screen. 


During the filling of the tank, vapor bubbles can be entrapped within the channels of 
the capillary device if the screen becomes wetted before the vapor can escape. Vapor can- 
not be permitted to remain within the channels of the device since it could cause drvout 
of the screen and failure of the ability of the device to expel gas-free liquid. Such va- 
por bubbles can be eliminated by pressurizing the tank, making the liouid subcooled with 
respect to the vapor pressure and causing the vapor to condense. Collapse of the vapor 
bubble must occur within a reasonable length of time (preferably minutes) so that the re- 
filling process can be completed and the subsequent mission for the spacecraft begun. 


A survey of the analytical and experimental investigations of bubble collapse can be 
found in reference 2. Of the work surveyed, that of Florschuetz and Chao3 seems the 
most ''oaprehensive. It defines the regimes in which inertia, heat transfer, or both mech- 
anisms determine the bubble collapse rate. For the case of heat transfer controlled col- 
lapse (of interest here) Florschuetz and Chao consider a solution based on the Plesset- 
Zwick temperature integral^ to be an uoper bound for the bubble sire versus time curve 
and their "plane interface" solution to be an "approximate lower limit". However, 
Prisnyakov® obtained a solution that "gives better agreement with experiments" and pre- 
dicts a faster rate of collapse. Likewise, the analysis of Theofanous^ predicts a fas- 
ter collapse rate than Florschuetz and '"■ao and "improvements in the agreement are noted" 
when non-equilibrium effects were consi :• red. 
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In applying any of the above theories to predict collapse tines durinq tank refilling 
there are two concerns. One is that the collapse tine for bubbles with volunes on the 
order of 100 cn^ is desired, while the above theories have only been verified with tests 
of bubbles on the order of 1 cn3. The second concern is that bubbles in contact with 
the inside walls of the channels will be elongated in shape, while the existing theory is 
applicable only to spherical bubbles. An analytical and experinental investioation was 
therefore perforned to deternine the influence of bubble size and shape on the collapse 
tine. 


Analysis 


Florschuetz and Chao3 define a dinensionless paraneter, B^ff, to classify the node 
of bubble collapse. For values of B^ff less than 0.05 heat transfer controls the col- 
lapse while values greater than 10 indicate inertia controls. An internediate case exists 
bettreen these values. The parameter is defined as 


B 


eff 



( 1 ) 


(A list of synbols can be found at the end of this paper.) 

Values of B^ff were calculated for typical conditions in a propellant tank and the range 


of conditions planned for the experinents, an 
is nuch less than O.OS (on the order of 10~^) 
trol the rate of bubble collapse. 


The configuraton of the bubble within the 
channel is as shown in Figure 1. Both 
screen and sheet metal surround the bubble 
on four sides with a vapor -liquid interface 
at each end. Under low-g conditions the 
liquid interface will have curvature, but a 
flat interface has been assumed here to sim- 
plify the analysis. It was assumed that the 
vapor bubble and liquid are initially in 
equilibrium and then the system pressure is 
instantaneously increased by some amount. 

The increase in pressure increases the sat- 
uration temperature of the vapor above the 
liquid temperature. This change in the 
thermodynamic condition results in condensa- 
tion of the vapor and collapse of the 
bubble. 


The change in the volume of the vapor 
bubble is dependent on the rate at which 
vapor condenses. 


dV - ^ dm (2) 


it was established that in all cases B^ff 
indicating that heat transfer will con- 



figure 1. Bubble Geometry 


The rate of condensation is dependent upon the rate at which heat is transferred i.om the 
vapor to the surrounding liquid. 


dm 



(3) 
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Therefore 


dv 




( 4 ) 


The 9 reater saturation tenperature causes the vapor to condense on the liquid, creating 
a liquid tile that is at the saturation tenperature. Condensation continues based on the 
rate at which this heat can be conducted into the liquid. Convection heat transfer is 
negligible since the vapor tenperature restains essentially unchanged during collapse. The 
unsteady heat conduction into a seni-inf inite solid is given by ^ 

Q - ^ (5) 


The cond> ctivity is that of the liquid only since the contribution of the thin sheet Mtal 
and screen has a neglibible effect on the heat transfer rate (this assunption will be dis~ 
cussed in sK>re detail later). Then 


iX. . kaxdt (6) 

A p L VWoT 
V • 


which can be reduced to 


A 





(7) 


where Ja, the Jacob number, is defined as 


Ja - 


ATCpP- 


LP., 


(8) 


Based on the bubble geometry in Figure 1, 


dV - ab dc 


(9) 


and 


A • 2ab + 2bc + 2ac 


(10) 


After integration, the following equation is obtained for the collapse of the bubble from 
its initial length (c^) to any final length (Cf ) . 


t 


4Ja^a 


ab 

2(a+b) 


In 


(a+b)Cj + abp 
(afb)Cj + abj 


( 11 ) 


For complete collapse ot the bubble Cf equals aero. It was assumed that only the length 
of the bubble changes as it collapses, but there would be a transition to a spherical bub- 
ble when the length c approached the channel thickness, a. Based on the assumptions it 
would be expected that this equation would be most applicable to the collapse of larger 
bubbles and be least accurate for the collapse of smaller bubbles and the final stages of 
collapse of any bubble. 
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It is interesting to note that if a spherical bubble geometry is used in solving equation 
(7 ) , then dV “ dr and the time for a bubble to completely collapse from an initial radius 
(rj^) is: 


t 


AJa^a 


(i:) 


This is the same result Florschultz and Chao3 obtained for their plane interface so)u- 
tion. Prisnyakov^ obtained a similar result except that the coefficient was 16 inst.iad 
of 4. 


For a given bubble volumet equations (11) ana (12) were used to calculate bubble col- 
lapse times. For larger bubbles the difference in geometry causes the rectangular shaped 
bubble to collapse about four times as fast as a spherical bubble of the same volume. 


Experiments 


In order to verify the analytical model presented in the previous section and to inves- 
tigate the influence of bubble geometry and the channel on bubble collapse, an experimen- 
tal investigation was performed. The approach was to form a bubble within a channel, 
pressurize the container in which the channel was installed and monitor the collapse of 
the bubble. Since a stationary bubble can be formed and confined within the liquid by the 
channel, a one-g test closely represents the low-g conditions. The mechanism of the bub- 
ble collapse, conduction heat transfer, is independent of the g-level and only minor chan- 
ges in the shape of the vapor bubble would be expected in low-g. 


A transparent channel, to permit viewing of the vapor bubble, was fabricated from plas- 
tic. It had an inside cross-section of 2.5 cm by 7.6 cm and was 3U.5 cm long. One side 
of the channel (see Figure 2) was a fine mesh screen having a 325 x 2300 (wires per inch 
in warp and shute directions) mesh, Dutch twill weave and an effective pore diameter of 7 
microns. This screen was capable of retaining any size vapor bubble within the channel. 
The channel was installed vertically within a transparent plastic box (Figure 3) . 


Freon 11 (CCI3F) was selected as the test liquid. This Freon has a boiling point of 
23.8°C at 1 atm so vapor bubbles could be easily created under ambient conditions. Liq- 
uid Freon 11 has the following proper ies at 20°C: p> 1.49 gm/cm3, l > 43.1 cal/gm, 

Cp « 0.205 cal/gm<K: ana o • 2.5 x 10~^m2/hr. The saturation curve is linear, 
having a slope of 0.32°K per kPa. 


The test procedure was to fill the channel and container with liquid so the channel was 
completely filled and submerged, excluding all air from the channel. While maintaining 
the container at a higher pressure to inhibit boiling, the channel was vented to form the 
vapor bubble. If necessary, a vapor generator could be used to aid in forming a bubble of 
the desired initial size, /he size of the bubble was monitored to ensure that the vapor 
was initially .n equilibrium with the liquid. The initial temperature of the liquid, sys- 
tem pressure and bubble length were recorded. 


A gaseous nitrogen pressurization system connected to the containet was set to give a 
desired increase in system pressure. The valve that applied the pressure increase to the 
container had an opening time that was negligible in comparison to the typical collapse 
times of 1 to 11 seconds. A motion picture camera photographed the bubble collapse. The 
oubble length versus time was measured using a scale on the channel and the frame rate of 
the camera. A total of 99 tests were performed, primarily varying the initial bubble 
length and the amount of pressure increase. 


For data correlation, the analytical model wa£ adapted to the specific test condi- 
tions. The plasti' walls of the channels influenced the modeling of the heat transfer, 
while the effect .ne screen could be neglected. The plastic was 1.3 cm 
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FiKure Mi'Ji‘1 


Kiguro i. Test Apparatus 


tniCK uno haa a ditfusivity ot 4.S x 10'^ m^/hr (twice the value for the liquid)) 
while the screen was O.Ub4 .-nm thick and liao a aiffusivity ot l.S x lO'^m^/hr (stain- 
less steel). The thermal (lenetration thickness given by' 


6 

t 




( 13 ) 


was used to evaluate the relative influence of the materials of the test article. For the 
vertical (elastic walls the heat conduction occurreo solely within the plastic so the liq- 
uid adjacent to the outside of the channel did not contribute to the heat transfer. The 
screen quicxly reached the temperature ot the surrounding liquid due to its thinness and 
high diffu..ivity, so its contribution to conduction perpendicular to its surface was neg- 
ligible. similarly, heat conduction parallel to the screen was shown to be negligible. 
Therefore, the heat transfer at the surfaces of the bubble were modeled as follows: 


upper surface and 3 vertical sides - unsteady conduction into plastic, and screen sur- 
face and liquid surface - unsteady conduction into liquid 


Equation (^) witn the appropr te values for k and A was used for both cases. 


From t.’ie film data it was established that the collapse of the bubble occurred as a 
rise 111 the liquid surface, changing only the length of the bubble, until small values of 
bubble length were reached. At bubble lengths less than one centimeter the bubble began 
tc dec;ease in width and during the last stages of collapse the bubble reached a 
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spherical shape and then disappeared. The data correlation concentrated on the initial 
sta 9 es of collapse when only the length of the bubble was chan 9 ing and the major change in 
volume occurred. 


As previously discussed it was established that the screen had a negligible effect on 
the heat transfer into the liquid. However, another potential influence of the screen is 
the effect of its flow resistance on the collapse rate. Liquid must flow through the 
.. screen, filling the channel as the vapor condenses. An analysis determined that the pres- 

% sure drop due to flow at the rate established by the bubble collapse had a negligible ef- 

fect on the pressure of the liquid within the channel. 


Preliminary correlations indicated that the bubble collapsed faster than predicted by 
equation (11) (including the above discussed modifications). Therefore a correlation co- 
efficient, F, was applied to equation (S) for the heat transfer rate, giving a term 
in the denominator of equation (11) . It was found that the value of F that best correl- 
ated the collapse time of the bubbles typically ranged from 1.3 to 2.0. Neither the ini- 
tial length of the bubble nor the change in system pressure appeared to have any effect on 
the variation in the value of F. Based on the excellent and more consistent correlation 
obtained for the bubbles having long collapi>e times, a value of 1.4 for F was selected as 
giving the best fit. This means that the coefficient in equation (11) is increased to 8, 
placing the value midway between that of Florschuetz and Chao3 with a coefficient of 4 
and Prisnyakov!^ with a coefficient of 16. 


Figure 4 is an example of the correlation of a test in which the bubble had a long col- 
lapse time. A very close match between the calculated and measured collapse rate was ob- 
taineo over the latter 10 seconds of the test. During the first 1.5 seconds of the test 
the bubble collapsed at a slower rate than predicted. This initial difference in the cal- 
culated and measured collapse rate becomes more evident in the shorter duration tests 
shown in Figures 5 and 6. Changing the heat transfer rate (through F) only changed the 
point at which the curves for the calculated ano measured collapse rate intersected and 
did not improve the match of their slopes. 




Figure <*, Correlation of Bubble Collapse, 
Long Collapse Time 


Figure 5. Correlation of Bubble Collapse 
Intermediate Collapse Time 
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Figure e. Correlation of Bubble Collapse, 
Short Collapse Time 



Figure 7. Collapse of Hydrogen Vapor Bubbles, 


Conclusions 


An analytical and experimental investigation of the collapse of a large vapor bubble 
inside of a channel has established that conduction heat transfer is the primary mechanism 
of collapse for the conditions of interest. It was found that the elongated shape of the 
bubble decreased the collapse time in comparison to a spherical bubble of the same vol- 
ume. An analytical model, based on conduction heat transfer to the surrounding liquid, 
gave excellent correlation of those tests having a longer collapse time ('vlO seconds). 


It appears that there are two stages to the bubble collapse, based on the tests per- 
formed here. In the initial stage, lasting about 1.5 seconds, the rate of collapse was 
less than predicted by the heat transfer model. Apparently the inertia of the liquid, 
flowing into the channel through the screen to replace the condensed vapor, reduced the 
initial collapse rate. This effect made the correlation of the shorter duration tests 
less accurate. 


During the later stage of collapse the rate was greater than predicted by the heat 
transfer model. This deficiency was corrected by the correlation coefficient which resul- 
ted in matching of the rates and times of collapse during that stage. This difference is 
typical of the variation noted in other bubble collapse analyses, as discussed in the in- 
troduction. When the collapse time is long, the influence of the initial phase became in- 
significant and the heat transfer model accurately predicted the collapse time. 

A correlation coefficient of 1.4 was selected as giving the best fit to all the data. 
Based on the data, the following equation will accurately predict the lonqer collapse 
times and it will predict too long a time for short collaps? periods. 


ab 

BJa^a 1^2 (a+b') (a+b)c^ + 


r 


(14) 
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When this equation is used to predict the collapse time of hydrogen vapor bubbles the 
result shown in Figure 7 is obtained. The saturation curve for hvdrooen is non-linear and 
the change in saturation temperature with pressure becomes small at pressures above 100 
kPa. For example, at 100 kPa a 50 kPa change In pressure causes the saturation tempera- 
ture to change by only l.S®*. Collapse times of many minutes, or even hours, are possi- 
ble with hydrogen. Hydrogen presents a "worst-case" in comparison to other propellants 
for the problem of collapsing entrapped vapor bubbles during tank refill. 


In comparison to the collapse time for a spherical bubble (based on Ref. 3), equation 
(14) yields a collapse time for an equal volume bubble in a channel that is about 10 times 
less. At these low collapse rates inertia effects should be negligible, so the assump- 
tions applicable to equation (14) are justified for this application and reasonably accur- 
ate predictions of the collapse time should be expected. 

List of Symbols 


A area 


a,b,c 

bubble dimensions 


®eff 

dimensionless group 


Cp 

specific heat of liquid 


F 

correlation coefficient 


Ja 

Jacob number 


k 

thermal conductivity of liquid 


L 

heat of vaporization 


m 

mass 


^P 

difference between system pressure 

and vapor pressure 

Q 

rate of heat transfer 


r 

bubble radius 


t 

time 


AT 

difference between vapor saturation 

temperature and 1 


temperature 
V volume 

o thermal diffusivity of liquid 

P liquid density 

Pv vapor density 

Pv average vapor density 

If temperature difference correction factor (see Ref. 3) 

2 thermal penetration thickness 
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